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^^ ' Abstract. We present GASPHOT, a tool for automated surface photometry and 

(<-^ \ morphological classification of galaxies in deep and wide fields. The requirements 

^vj . for any such tool are reviewed, and its use for the classification of high-z galaxies is 

' presented. In the case if HDF-like images, for galaxies having a magnitude ranging 

jr~^, from 24 to 27.5, the uncertainties on the photometric parameters derived from 

^ . GASPHOT are respectively AAf = 0.02 - 0.1, Alogi?^ ~ 0.03, An = 0.02 - 0.5. 

fT^ ' Keywords: methods; data analysis — galaxies: photometry — galaxies: classifica- 

, tion. 
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ff^ , 1. Introduction 

O; 

Q ■ In the last years the Hubble Space Telescope provided very deep images 

^^ , of field galaxies, giving a glimpse of their very first evolutionary phases. 

^^ I We now have a chance to understand the basic elements of galaxy 

,J^ ■ formation and evolution, but we need new instruments to handle the 

^p^l thousands of very faint objects that could be present in one single 

O ' image. 

In particular, it has been noticed that at high redshift - and for this 



U 



^ \ kind of images - the Hubble classification system could be, if not totally 

wrong, at least not practical to use (Abraham et al., 1996). On the one 
hand, the objects in the images are small, with a typical radius of a few 
^ ' pixels. Thus, it is very hard to detect the fine structure elements needed 

j^ ■ to distinguish the different classes of spirals or to separate the barred 

and unbarred families of galaxies (but see Abraham et al., 1999). On the 
other hand, the morphological statistical analysis of the Hubble Deep 
Fields (Abraham et al., 1996) showed that at high redshift (z^O.T) the 
percentage of Peculiar and Irregular galaxies could be very high (up to 
30%). 

Thus, the Hubble diagram does not appears very useful to distin- 
guish between the different classes of high redshift galaxies. On the one 
hand, it brings plenty of details for distinguishing between sub-classes 
of galaxies (the spirals) that can not be separated for distant objects; 
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on the other hand, it gives few informations about a class of galaxies 
(the Irregular) which is dominant at high redshift, and that we would 
like to split at least in truly irregular, peculiar and interacting galaxies. 

For these reasons, different authors tried to build new classification 
systems, based on different quantitative parameters. The most used is 
the luminosity profile, measured either with a "concentration param- 
eter" (Abraham et al., 1994) or with a slope of the surface brightness 
profile; this can also be combined with colors, asymmetry (Abraham 
et al., 1996; Conselice et al., 2000) or power spectrum in Fourier space 
(Takamiya, 1999). 

In order to extract the valuable photometric parameters (total mag- 
nitude, optical radius, b/a, morphological type T) for the thousands 
of objects expected in large fields, we can not relieve on the usual 
photometric tool, which are optimized to be used on single objects. Such 
tools always involve some amount of interactivity, which is ruled out 
by the large number of galaxies now present. We need the photometric 
analysis process to be completely automatic. 

While there are many tools available for the detailed surface pho- 
tometry of single galaxies, there are few instruments which are build 
for the study of large databases of objects (and they are not for public 
use). Usually the extraction and "aperture photometry" is handled by 
tools such as FOCAS (Jarvis and Tyson, 1981) MORPHO (Odewahn, 
1995) or SExtractor (Bertin and Arnouts, 1996), while the only soft- 
ware developed for automated surface photometry of a large number 
of galaxies are GIM2D (Marleau and Simard, 1998) and the HST MDS 
software (Ratnatunga et al., 1999). 

2. Structure and performances of GASPHOT 

For the reasons described above, we started creating a photometric tool 
for detailed surface photometry of large images. The process is mainly 
divided in two steps: first, a modified version of SExtractor (Bertin 
and Arnouts, 1996) is used to perform the identification and aperture 
photometry of each object. Then a second program is used to analyze 
the photometric profiles and derive the main photometric parameters 
of each object. Taking into account the effects of the convolution with 
the PSF, the aperture photometry is fitted with a Sersic law fi oc r^'^ 
with five free parameters: the total magnitude Mtot, the half-luminosity 
radius Re, the Sersic index n (which will be taken as a morphological 
index), the flattening b/a and the value of the local background. 

The bias and errors in the determination of parameters have been 
estimated by running IRAF simulations of galaxies adopting the read- 
out noise, gain, background level and PSF of the HST WFPC2 detector. 
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Morphology of high-z galaxies o 

We performed simulations with galaxies having magnitude from 23 
to 27.5, and optical radii from 3 to 10 pixels. First, in order to re- 
move the problems due to blending effects, we analyzed the results of 
GASPHOT on a sample of about 2000 galaxies positioned over a grid. 
Galaxies were assumed to follow a pure exponential or de Vaucouleurs 
law, but we also made a few tests with galaxies following a Sersic 
law with n ranging up to 6. For galaxies having a magnitude of 24- 
27, we obtained uncertainties of AM = 0.02 - 0.1, Alogi?e = 0.03, 
An = 0.02 — 0.5, with no sizable bias. Using the Hubble Deep Field im- 
age parameters, the limit magnitude to obtain meaningful informations 
from the photometry seems to be around 27.5. The computational time 
seems acceptable for images with thousands of objects (about 20-30 
seconds/object). 

The morphological classification of the objects has been done on the 
basis of the Sersic index n that best fits the observations. In the future, 
we plan to include an asymmetry parameter. While the separations 
between elliptical and spiral galaxies seems to be excellent (see Fig.l) 
we still have to test the behaviour of the software when dealing with 
mixed types such as SOs. 

As a second, harder test, we spread 500 galaxies uniformly with a 
random distribution over a 1400x1400 pixel image. The test included 
many galaxies heavily blended or with a close companion. We also 
pushed our magnitude limit up to 28.5, and tested our tool to reproduce 
the original parameter of the galaxies. For comparison, we tested the 
most common software, SExtractor, on the same image. We only could 




Figure 1. Morphological classification of galaxies on grid, for galaxies having total B 
magnitude of 27 in B. Note that early- and late- type galaxies are clearly separated. 
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Figure 2. Results for a simulated image of 500, uniformly distributed, elliptical 
and spiral galaxies. On the left side we show the disagreement of the magnitudes 
measured by SExtractor with the "real" magnitudes of the simulated galaxies. Full 
dots represent elliptical galaxies, empty squares spiral galaxies. The same test is 
performed on the right, for the same image, by GASPHOT. 



perform the test for magnitudes in this last case, because SExtractor 
does not provide a morphofogical classification or the radius Re (Fig. 2) . 
In the future, we plan to produce a morphological and photometric 
catalog of the Hubble Deep Fields using GASPHOT, and to extend 
its application to a sample of low-to-intermediate redshift clusters. We 
also plan to produce a public release of the software to be available in 
the next few months in our web page. 
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